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ABSTRACT 



£ — ' Context. The impact of rapid rotation on stellar evolution theory remains poorly understood as of today. Vega is a special object in 

£Nj , this context as spectroscopic and interferometric studies have shown that it is a rapid rotator seen nearly pole one, a rare orientation 

particularly interesting for seismic studies. In this paper we present a first systematic search for pulsations in Vega . 
Aims. The goal of the present work is to detect for the first time pulsations in a rapidly rotating star seen nearly pole-on. 
Methods. Vega was monitored in quasi-continuous high-resolution echelle spectroscopy. A total of 4478 spectra were obtained within 
, 3 individual runs in 2008, 2009 and 2010. The resolution was higher than R = 65000. This data set should represent the most extensive 

high S/N, high resolution quasi-continuous survey obtained on Vega as of today. Equivalent photospheric absorption profiles were 
calculated for the stellar spectrum, but also the telluric lines present in the spectrum, the latter were used as a radial velocity reference. 
»-^H' Residual velocities were analysed and periodic low amplitude variations, potentially indicative of stellar pulsations, detected. In 

a subsequent step, LSD profile variations were searched for in a bidimensional analysis night by night. Finally, a Lomb-Scargle 
j.^ ■ periodogram of each velocity bin of the profile was performed for all three observing runs. 

Results. Based on high resolution echelle spectroscopy, we have obtained indications of periodic variations of very small amplitudes 
within the residual radial velocity curves of Vega. All three data sets revealed the presence of residual periodic variations: 5.32 and 
9.19d~' (A » 6ms- 1 ) in 2008, 12.71 and 13.25 d" 1 (A » 8ms-') in 2009 and 5.42 and 10.82 d^ 1 (A « 3-4ms-') in 2010. However, 
it is too early to conclude that the variations are due to stellar pulsations, and a confirmation of the detection with a highly stable 
spectrograph is a necessary next step. 
^ ' Conclusions. If pulsations are confirmed, their very small amplitudes show that the star would belong to a category of very "quiet" 

pulsators. 

O ■ 

^ . Key words, stars: oscillations - stars: rotation - stars: individual: Vega - asteroseismology 

■ 1. Introduction field no tably in the description of the organization of p-modes 

spectra (Ligni eres & Ge oreeot, 2009). Even if we understand the 

Rapid rotation is a characteristic feature of massive and interme- basic properties of the frequency spec trum, the mode identifica- 

diate mass stars that has a significant and yet poorly understood tion process remains a difficuk task In particular> in the p . mo de 

impact on their structure and evolution. It is a current challenge frequency ^ the large rotationa i sp i it tings between modes 

^ ; of stellar physics to obtain observational constraints that will im- of difkrmt azimuth al number mix the frequencies of different 

. prove the modeling of the physical processes induced by rota- multiplets . Qne possible solution to this problem would be to ob- 

■ tion. Seismology should play a particular role in this context as serye a rapidly rotating pulsating star seen pole . on Cancellation 
it is the only way to directly probe stellar interiors. effect - m d i sk .i ntegrat ed light would then select axisymmetric 

At present, however, the oscillation spectra of rapidly rotat- mode ^ non-axisymmetric modes being cancelled out. Such a 
ing stars can not be unambiguously interpreted. One difficulty election effect would considerably simplify the observed spec- 
resides in the modeling of the effects of rapid rotation on the trum and thus the mode identification process. 
osciUation modes. There has been some recent progress in this Vega ifj knQwn tQ fee a rapidly rotating star and reveak a 

; nearly pole on inclination. This has been shown by analyzing the 

Send offprint requests to. T. Bohm spectroscopic and interferometric signatures of the surface grav- 

* Based on observations obtained at the Bernard Lyot Telescope h darkening ( Gulliver et all [T9MlPeterson"eTaIl |2006). The 
(TBL, Pic du Midi, France) of the Midi-Pyrenees Observatory, which , , • . c . • , • , • •, 

' , , , t ■ xt ■ i j o • j „„ • /tx to t t\ spectroscopic and interferometric approaches provide a similar 

is operated by the mstitut National des Sciences de 1 Umvers (1NSU) . .. . , . , _ , rf , , . , 

r .£ „ . XT t - i A i r> u u c ■ tc /Ax.™ <^ inclination angle l = 5-7 degrees, although they diner in the re- 

ot the Centre National de la Recherche Scientifique ot France (CNRS), b , . b , , _i 

and at the Canada-France-Hawaii Telescope (CFHT), which is oper- P orted equatorial velocities, v eq = 245 + 15 kms andv e q ~175 
ated by the National Research Council of Canada, INSU/CNRS and the km s~' from the s pectro scopic analysis of IGulliver et alj (Il994h 
University of Hawaii. and iTakeda et al.l (2008), respectively, and a significantly larger 
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value , v ea ~ 275 kms ' , from interferometry (iPeterson et all 
2006; Aufdenberg etafl 120061) . As a seismic target, Vega also 
presents the important advantage that its fundamental parame- 
ters are very accurately determined, particularly when compared 
with other rapidly rotating stars. 

As reviewed bv lGravl(l2007h . Vega was used since more than 
150 years as a photometric and spectrophotometric standard star, 
the historic reason being most likely its brightness and near- 
zenith position for european and north american observatories. 
This particular role of Vega can be considered surprising, since 
several authors report photometric variability of at least 1 - 2%, 
but with no evidence of periodic variations. A previous study of 
spectr oscopic variability has been conducted bv lGorava & S ingh 
( fl983h who described a variabil ity of Ha alfhought it w as not 
confirmed in subsequent studies (Charl ton & Meverlll985h . 

The location of Vega in the HR diagram is roughly that of 
main-sequence AO stars, but bluewards of the 6 Scuti instability 
strip, ther efore stellar os cillations are in principle not expected. 
However, iFernie I d!98lh concluded that historic radial velocity 
curves and photometric variability tend to indicate episodical 5 
Scuti type pulsation, with a frequency close to the fundamental 
radial order of Ff un d ~10d~' for Vega's parameters as known at 
that time. The most recent determ ination of fundamental stellar 
parameters bv lTakeda etal] d2008l) led him to conclude that this 
rapid rotator has parameters strongly varying from pole to equa- 
tor, more precisely an effective temperature from T e g = 9867 to 
8931 K, a surface gravity from \ogg = 3.308 to 3.313 (where g 
is expressed in cm s~ 2 ) and a radius from R = 2.520 to 2.763 R Q , 
respectively. 

Until know, no rapidly rotating, pole-on and observably pul- 
sating star was known. The goal of our work was therefore to 
search for low amplitude stellar pulsations in an extensive high- 
resolution quasi-continuous spectroscopic survey of Vega. 

Sect.|2]describes the observations and data reduction, Sect. [3] 
presents the analysis of radial velocity variations in the spectro- 
scopic time series of Vega, and in Sect. [4] results are discussed 
and a conclusion is drawn. 



2. Observations and data reduction 

The analysis presented in this paper is based on three datasets 
obtained in 2008, 2009 and 2010. Vega was monitored dur- 
ing 3 nights in 2008 (26 th , 27 th and 29 th of July) with the 
high-resolution echelle spectropolarimeter NARVAL at the 2m 
Bernard Lyot Telescope (TBL, Pic du Midi, France). The first 
two nights observations were acquired at R = 65000 in polari- 
metric mode, and the last night at R = 75000 in spectroscopic 
mode. We obtained 1213 spectra covering 19.9 hrs. In 2009 we 
obtained 1293 spectra covering 13.7 hrs on 3 nights (9 th - 11 th 
of September 2009) at ESPaDOnS/CFHT (R = 68000) in polari- 
metric mode. Finally, in 2010, we got a huge data set of 1972 
spectra covering 28.8 hrs on 5 nights (15 th - 19 th of July 2010), 
obtained again in the polarimetric mode at NARVAL/TBL. Table 
Q]summarizes the log of the observations during these 3 runs. 

The general observing strategy was to obtain as many maxi- 
mum S/N observations of the target during the night as possible. 
The goal was to obtain a versatile data set, enabling us to search 
for the presence of a magnetic field as well as stellar pulsations. 
Results of the s pectropolarimetr i c stud y of the first observing run 
are presented in Lignier es et alJ J2009), where the detection of a 
weak magnetic field in Vega has been announced, revealing the 
existence of an as-yet unexplored class of weakly magnetic stars 
in the intermediate mass domain. The joint analysis of the two 



first runs (2008 and 2009) confirmed this first detection and indi- 
cated a rotational modulation of the polarization signature with 
a period of 0.732 ± 0.008 d, a value which is very close to the 
rotation period announced by Takedaet al. (2008]) (^ 0.733 d). 
These results are published in Peti t et alJd2010l) . A recent analy- 
sis of the 2010 d ata set yielded a sli ghtly different rotation period 
of 0.678 + g : g^ d dAlina et all 1201 lh . 

The global data set containing 4478 spectra obtained at the 
highest possible signal to noise should represent, as of today, the 
most extensive high resolution quasi-continuous spectroscopic 
monitoring of Vega. It should be noted that for such a bright 
star the exposure times are extremely short (4 to 13 seconds 
on the star) in order to avoid the saturation of the CCD, while 
the readout times of the CCD dominate the duty-cycle, the aver- 
age time interval between two exposures being close to 0.8 min. 
Subexposures obtained in polarimetric mode provide as a direct 
result from the data reduction pipeline, in addition to the polar- 
ization information, a perfectly calibrated intensity spectrum. 

The wavelength covered extends from 369 to 1048 nm for 
both spectrographs. Standard calibration spectra were obtained 
at the beginning and at the end of the night; most importantly 
the absolute wavelength calibration relies on a Th/Ar spectrum 
at the beginning of the night. Most of the data reduction was car- 
ried out following standard reduction procedures using a dedi- 
cated spectroscopic data reduc t ion pa ckage, based on the algo- 
rithm described in lDonati et all dl997l) . T his package al so makes 
use of the "optimal extraction algorithm" dHornel[l986f) . A stan- 
dard heliocentric velocity correction was performed. The in- 
trinsic wavelength calibration accuracy achieved with the 2D- 
polynomial fit procedure is better than 150 ms -1 for both data 
sets. In order to be able to calibrate potential spectral shifts dur- 
ing the night a cross-correlation of telluric lines present in the 
spectra was performed by the standard reduction package. 

After all shifts have been applied to the wavelength grid, 
the residual telluric line velocity should reflect precisely the op- 
posite of the heliocentric velocity correction. In order to ver- 
ify this, we comput ed least square deconvolved (LSD) profiles 
dDonati et all fl 997) of the more than 100 telluric water vapor 
lines contained in each stellar spectrum, and determined their 
centroid by fitting a gaussian profile. We discovered the presence 
of residual radial velocity errors of the telluric lines. In order 
to optimize the radial velocity precision we corrected then for 
the remaining velocity shifts. The residual spectral stability is 
around 10-15 ms -1 (depending on the night and on the site), tak- 
ing into account all the multiplex information of the lines present 
in the spectra. 

At this point it should be noted that a stability ana l ysis o f 
O2 telluric lines has been carried out by ICaccin et alJ (Q985). 
They report that the line positions change more or less errati- 
cally with time, their displacement depending upon the horizon- 
tal wind structure above the observing site. They note however, 
that unless the observations are made near the horizon, the com- 
ponent of the wind velocity along the line of sight is very seldom 
significant, the displacement reaching up to 10 ms~' only in very 
severe conditions (large air mass, strong winds). In low wind 
conditions and for only one to three airmasses their time insta- 
bility should be probably of the order of 5 ms~' . This is the case 
for our observations of Vega, since they have been obtained in 
airmasses between 1 and 2, and only during the last night of the 
CFHT run airmasses up to 2.6 were reached. In a recent work, 
Figu eira et alJ (l2010h evaluated the stability of atmospheric lines 
with HARPS (on the 3.6 m telescope at LaSilla, ESO), and con- 
cluded in a similar way to lCaccin et alJ d 19851) that a precision of 
5 ms -1 can routinely be obtained in restricted airmass conditions. 
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Table 1. Log of the spectroscopic (& polarimetric) observations of Vega. (1) Date of the observation (UT) dPetit et al.l (1201 Oh 
indicates the date in local time, leading to a -1 day shift for the Sept. 2009 CFHT data). (2) and (3) Heliocentric Julian date (mean 
observation, 2,450,000+) of the first and the last stellar spectrum of the night, respectively; (4) Number of high resolution RS Cha 
spectra; (5) exposure time (sec);(6) total hours covered on the sky; (7) nightly average and standard deviation of signal to noise per 
resolved element at 520 nm. 



Date 


JDfirst 


JD last 


N 

1 " spec 


texp (SCC) 


t cov (hrs) 


S/N 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Jul. 26 2008 


4674.3489 


4674.6072 


389 


6 


6.2 


812+50 


Jul. 27 2008 


4675.3748 


4675.6551 


390 


6 


6.7 


822±61 


Jul. 29 2008 


4677.3480 


4677.6415 


434 


6 


7.0 


823+82 


Sep. 9 2009 


5083.7400 


5083.9066 


402 


4 


4.0 


865+85 


Sep. 10 2009 


5084.7111 


5084.8787 


331 


4 


4.0 


887±69 


Sep. 11 2009 


5085.7122 


5085.9498 


560 


4 


5.7 


736+155 


Jul. 15 2010 


5393.3643 


5393.5735 


355 


13 


5.0 


1117±37 


Jul. 16 2010 


5394.3635 


5394.6576 


496 


13 


7.1 


1051+158 


Jul. 17 2010 


5395.3650 


5395.6510 


453 


13 


6.9 


1197+59 


Jul. 18 2010 


5396.3580 


5396.4662 


185 


13 


2.6 


1222±16 


Jul. 19 2010 


5397.3561 


5397.6569 


483 


13 


7.2 


1170±70 



What we retain from this latter work is the fact that increasing 
airmass induces asymmetry of telluric lines, this asymmetry is 
therefore expected to have an impact on the measured RV. This 
behaviour-depending on the airmass of the observed object re- 
peats therefore every night for the same target (within a 3 or 5 
night observing run). Any periodic behaviour with different peri- 
odicities than Id cannot be attr ibuted to this i ntrin sic telluric line 
variability. Moreover, even if ICaccin et al .1 (Il985h invoke high 
altitude atmospheric movements or pressure shifts as responsi- 
ble for inducing asymmetry in telluric lines, we could not find 
any report of sinusoidal radial velocity variation of telluric lines 
during an observing night. 



The next step of the data reduction was to calculate for all 
stellar spectra photospheric LSD profiles, using a mask corr e- 
sponding to the spectral type AO (see also jPetit et al. | (|20 10)). 
We ha ve to notice that, as lGulliver et al.l J 19941) and Takeda et al.1 
(2008) have shown, the spectrum of this rapidly rotating star 
is very inhomogeneous, since the weak spectral lines form ex- 
clusively in the almost 1000K colder gravity darkened equato- 
rial zone, leading to a flat-bottomed trapezoidal profile shape. 
However, the more numerous strong lines, and consequently the 
mean LSD profile, come from intermediate longitude to polar re- 
gions; therefore an AO mask represents a good choice for Vega. 



The equivalent photospheric profiles contains multiplex in- 
formation from more than 1100 individual photospheric lines. 
Typical S/N values of the resulting LSD intensity (Stokes I) pro- 
files are consistently around 1000 per resolved eleme nt, follow- 
ing the behaviour pointed out by Don ati et~ai1 d 1997b . and indi- 
cating that the convolution model underlying the LSD cannot be 
trusted above this level of accuracy (this limitation is not valid 
for the other St okes profiles) . Desp ite its almost pole-on geom- 
etry of i = 7° dTakeda et all 120081) . the LSD profiles of Vega 
(vsin/ = 22 kms" 1 ) needed to be fitted by rotationally broad- 
ened profiles in order to perform an accurate radial velocity de- 
termination. The finaly achieved precision in radial velocity is 
again 10-15 ms -1 for each stellar LSD profile of this data set. 
The only way to take into account a potential instrumental drift 
during the night is to work on the difference between stellar and 
telluric lines, knowing the limitations of the telluric lines as an 
absolute radial velocity reference. 



3. Search for periodic variations in Vega 

3.1. Radial velocity variations in Vega 

Due to the one year time gap between the three runs we de- 
cided to analyse each run separately. Plotting each radial velocity 
time series separately, we observed the presence of trends within 
some nights of less than 100 ms^'/mglit; for the TBL 2008 run, 
and less than 60 ms 1 /night for the CFHT 2009 and TBL 2010 
run. These slow variations did not appear to be sinusoidal, and 
we decided to take out these largescale trends by normalising the 
radial velocity curve of each night by a 2nd degree polynomial. 
In order to check for consistency, we analysed all three couples 
of normalized and raw radial velocity curves and concluded on 
the coherence of this approach by detecting the same frequencies 
in the range above 5 d~' . 

Figs.[T]|2]and|3]show the residual radial velocity variations of 
Vega observed with NARVAL/TBL (2008), ESPaDOnS/CFHT 
(2009) and NARVAL/TBL (2010), respectively. These curves 
represent the radial velocity variations after the nightly polyno- 
mial normalisation. Superimposed are the results of the corre- 
spon ding frequency analysi s (Table |2]i with the Period04 pack- 
age dLenz & Breger L l2005h . This package is based on the dis- 
crete Fourier Transform algorithm and multi-sine least-square 
fits of the data. Interestingly, night 3 of the TBL (2008) data set 
was the only night in which data were acquired in pure spec- 
troscopic mode; Fig.[T]shows well that the radial velocity curve 
does in no way depend on the acquisition mode, since variations 
are of the same order of magnitude and in phase within all three 
nights. This rules out the polarimetric mode as generator of faint 
radial velocity variations. 

Confiden c e lev els have been calculated according to 
Breger et al. ( 1993), who studied empirically the possibility of 
a peak in the periodogram being a true signal of pulsations 
with respect to t he noise level; this work was refined later by 
lKuschnigl(ll997l) : the significance of a peak in the amplitude pe- 
riodogram exceeding 4.0 times the mean noise amplitude level 
after prewhitening of all local frequencies has a 99.9% proba- 
bility to be due to stellar pulsations (99.0% for a ratio of 3.6, 
90.0% for a ratio of 3.2). Figs, g] [5] and [6] show the results of 
the Period04 frequency analysis for the three different runs. The 
two dominant frequencies in the data set of 2008 are FI2008 = 
9.19 d 1 (A = 5.8ms- 1 )andF2 20 08 = 5.32d" 1 (A = 5.9 ms- 1 ). 
We stopped the frequency analysis after extraction of the two 
strongest frequencies, the next potential frequency having a sig- 
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nificantly fainter amplitude. In the 2009 data set, the two dom- 
inant frequencies are FI2009 = 12.71 d _I (A = 8.2 ms _1 ) and 
F22009 = 13.25 d -1 (A = 7.8 ms~'), the next frequency shows an 
amplitude of around 3.6 ms _1 and was not taken into account. 
Finally, in the 2010 data set variations are even fainter in ampli- 
tude with FI2010 = 5.42 d -1 (A = 4.2 ms _1 ). All frequencies of 
the different data sets are summarised in Tab. |2j we added a po- 
tential second very low amplitude frequency from the 2010 data 
set for comparison purposes. 

TBL 2008 observations 




4674.35 4674.40 4674.45 4674.50 4674.55 4674.60 4674.65 




4675.35 4675.40 4675.45 4675.50 4675.55 4675.60 4675.65 




Fig. 1. Residual radial velocity variations of Vega 
(NARVAL/TBL) 2008. Superimposed (continuous line) is 
the result of the corresponding frequency analysis (Table 13. 
Time is expressed in HJD = 2450000 + A HJD. 



CFHT 2009 observations 



50 

40 - 
20 - 
- 
-20 - 
-40 - 
-60 L 




o o 

5083.75 5083.80 



5083.85 5083.90 5083.95 5084.00 




5084.70 5084.75 5084.80 5084.85 5084.90 5084.95 5085.00 




5085.70 5085.75 



5085.80 5085.85 

AHJD 



5085.90 5085.95 5086.00 



Fig. 2. Residual radial velocity variations of Vega 
(ESPaDOnS/CFHT) 2009. Superimposed (continuous line) 
is the result of the corresponding frequency analysis (Table 
Time is expressed in HJD = 2450000 + A HJD. 



TBb 2010 observations 



60 
40 
20 


-20 
-40 




5393.35 5393.40 5393.45 5393.50 5393.55 5393.60 5393.65 




5395.35 5395.40 5395.45 5395.50 5395.55 5395.60 5395.65 



60 
40 
20 


-20 
-40 
-60 



5396.35 5396.40 5396.45 5396.50 5396.55 5396.60 5396.65 




Fig. 3. Residual radial velocity variations of Vega 
(NARVAL/TBL) 2010. Superimposed (continuous line) is 
the result of the corresponding frequency analysis (Table [2}. 
Time is expressed in HJD = 2450000 + AHJD. 



Vega Fourier Periodogram TBL 2008 




20 30 

Frequency (c/d) 



Fig. 4. Result of a Period04 frequency analysis on the 
NARVAL/TBL data of 2008. Confidence levels are: 99.9% 
(cont), 99.0% (dotted) and 90.0% level (dashed). These lev- 
els were level calculated on a 10 d 1 sliding window after 
prewhitening with Fl and F2. 



3.2. Line profile variations in Vega 

For almost pole-on objects like Vega, radial velocities are par- 
ticularly sensitive to zonal modes of low degree £, if expressed 
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Vega Fourier Periodogram CFHT 2009 




10 20 30 

Frequency (c/d) 

Fig. 5. Result of a Period04 frequency analysis on the 
ESPaDOnS/CFHT data of 2009. Confidence levels are: 99.9% 
(cont), 99.0% (dotted) and 90.0% level (dashed). These lev- 
els were level calculated on a 10d~' sliding window after 
prewhitening with Fl and F2. 



Vega Fourier Periodogram TBL 2010 

n = 5.42 c/d, A = 4.2 m/s 




20 30 
Frequency (c/d) 

Fig. 6. Result of a Period04 frequency analysis on the 
NARVAL/TBL data of 2010. Confidence levels are: 99.9% 
(cont), 99.0% (dotted) and 90.0% level (dashed). These lev- 
els were level calculated on a 10 d~' sliding window after 
prewhitening with Fl. 

Table 2. Frequencies and amplitudes of the main frequencies 



ID 


freq. 


A 




d" 1 


m s -1 


FI2008 


9.19 


5.8 


F22008 


5.32 


5.9 


F12009 


12.71 


8.2 


F22009 


13.25 


7.8 


F120I0 


5.42 


4.2 


F2 2 oio 


10.82 


3.3 



in traditional spherical harmonics. We wanted to see if any 
periodic variation corresponding to a mode of higher degree 
I could be seen in the equivalent photospheric profiles them- 
selves. We therefore analysed the temporal line profile variations 
of the longest data set (2010) for each individual velocity bin 
of 1.8 km s by comp uting a 2D (velocity-frequency) Lomb - 
Scargle periodogram (Scargle, 1982; Home & Baliunas 

, 1986), 

shown in Fig. [7] after having subtracted the average profile of the 
run. It should be noted that line profile variations during this run 
are fainter than AF/F C < 7x10 4 , where AF and F c are the profile 
variations and the (normalized) continuum flux, respectively. 



No indication of the presence of any particular frequency 
within the profile, which was not also present in the continuum, 
is detected. The weaker signal of the 2D periodogram within the 
line profile is due to the square root dependency of the noise, 
diminishing the absolute noise amplitude within the absorption 
profile. Fig.|8]shows the quadratic sum of the individual Scargle 
periodograms of all velocity bins within the continuum versus 
the sum of those within the line. Despite the logarithmic repre- 
sentation, no noticeable difference in the line- versus continuum- 
Scargle periodogram can be seen. 



Scargle 2D Periodogram Vega TBL 2010 




velocity (km/s) 



Fig. 7. Two-dimensional Scargle line profile analysis, corre- 
sponding to a frequency search for each velocity bin. Data set 
of 2010 TBL (5 nights). The residual line profile variations have 
a standard deviation of 0.0007 F/Fc. Dashed lines indicate the 
+v sin /limitations as well as the centroid of the photospheric 
equivalent profile. 



4. Discussion and conclusions 

We searched for pulsations in Vega as a near pole-on config- 
uration is particularly favorable for the seismology of rapidly 
rotating stars. 

Equivalent photospheric line profiles were calculated, using 
the multiplex information of all photospheric spectral lines in 
the observed wavelength range. In this instrumental configura- 
tion, the highest attainable radial velocity calibration was ob- 
tained by taking the telluric water vapor line profiles present in 
the spectrum as absolute wavelength references. A subsequent 
radial velocity analysis with Period04 was performed and low 
amplitude (A < 8m s~') multi -periodic variations in the residual 
radial velocity curves were detected, with frequencies ranging 
roughly from 5 - 13d 1 , depending on the run. The frequencies 
detected in this work, if confirmed, are consistent with values 
occasionally r eported in the past, as described in the review by 
iFernie I (119811) . 

The observed amplitudes of the variations are very small 
with respect to the intrinsic calibration precision, but taking the 
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Scargle LSD — profile periodogram TBL 2010 
1 00-00 F I 1 1 I I J| 1 1 I 




10 20 30 40 50 

frequency (c/d) 

Fig. 8. Equivalent Scargle periodogram for the line (continous) 
and the continuum (dotted) during the TBL 2010 run, showing 
that line and continuum have the same frequency pattern, and 
therefore no indication of periodic line profile variations are de- 
tected at this level of precision. 



telluric water vapor lines of the spectrum as a reference en- 
abled us to reach a precision on the individual radial velocity 
measurement of 10-15 ms -1 . This auto-calibration of the spec- 
trum should rule out, in principle, any instrumental origin. While 
the spectroscopic and polarimetric set-up of NARVAL/TBL pro- 
duces residual radial velocity variations of the same amplitude, 
and in phase with each other, the configuration of the spec- 
trograph itself can quite certainly be ruled out as an origin. 
Nevertheless, the similarity of frequencies observed in both runs 
at TBL (2008 and 2010), quite different to the frequencies ob- 
served in 2009 at CFHT, leave us somewhat suspicious concern- 
ing a potential, not yet understood, instrumental origin. 

The small amplitudes of the detected periodic variations 
were obtained after detrending the individual nights of each run 
with 2nd order polynomials. We must keep in mind that these 
nightly drifts of the radial velocity were less than 100ms _1 in 
the 2008 run, and less than 60 ms 1 in the CFHT 2009 and TBL 
2010 run, and a priori non periodic. No origin could be at- 
tributed to these variations. The positive point however was that 
we were able to detect the periodic signals in both, the raw and 
the nightly-de trended radial velocity curves. 

The intrinsic telluric line variations due to strong, high al- 
titude wind conditions and high airmasses can be of the same 
order of magnitude as the detected periodicities. These extreme 
values of intrinsic shifts were certainly not reached, since air- 
masses during our observations were always small. Moreover, 
and as an important point, the coherence of the detected periodic 
signals during 3 or 5 observing nights (depending on the run) 
and the fact that the variations are not in phase with the begin- 
ning of each night tend to rule out the telluric lines as origin of 
the periodic variations, and strengthen the interpretation as small 
amplitude stellar pulsations. 

The detected frequencies are of the same order as the fun- 
damental mode calculated for Vega. Nevertheless, an excitation 
mechanism for such modes in this region of the HR Diagram is 
not known yet. 

At this stage we can ascertain that should Vega be confirmed 
as a pulsator in necessary subsequent high-precision observa- 



tions with highly stable spectrographs, such as Sophie/OHP, the 
very small amplitudes observed during the three runs in 2008, 
2009 and 2010 show that the star would belon g to a categor y 
of very "quiet" pulsators. Still, as suggested by Fernie (1981), 
some rarely occuring, episodic higher amplitude pulsations can 
not be ruled out. If pulsations were confirmed in the future, Vega 
would certainly be one of the most interesting stars to monitor 
in high-precision spectroscopy and photometry. 

Another approach to search for rapidly rotating, pulsating 
and pole-on orientated stars is to search for such candidates 
within the (5-Scuti instability strip; this work has been recently 
initiated by our group. 

As a conclusion, one thing is certain: Vega has revealed these 
past years to be a more and more enigmatic star and clearly de- 
serves our further attention! 
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